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ABSTRACT 


Electromagnetic  scattering  from  rough  surfaces  is  of  prime 
importance  in  the  engineering  field  since  it  affects 
communications,  radar,  remote  sensing,  acoustics,  etc.  The 
actual  problem  of  scattering  from  rough  surfaces  is 
complicated  and  involves  three  dimensional  scattering  from 
either  lossy  or  dielectric,  electrically  large  surfaces. 
Integral  equations  are  widely  utilized  to  solve  this  kind  of 
problem  but  this  solution  to  the  problem  is  generally 
computationally  intensive. 

In  the  On-Surface  Radiation  Boundary  Condition  (OSRC) 
method,  a  higher  order  radiation  condition  is  imposed  directly 
on  the  surface  of  the  scatterer.  This  reduces  the  integral 
equation  for  the  scattered  field  to  a  line  integral  which  can 
be  easily  evaluated  numerically. 

In  this  thesis,  the  OSRC  method  is  used  to  formulate  the 
problem  of  scattering  from  periodic  rough,  two-dimensional 
surfaces  illuminated  by  a  transverse  magnetic,  p''ane 
electromagnetic  wave.  Three  geometric  surfaces  are 
considered.  A  comparison  is  made  between  tho  present 
formulation  the  exact  solution,  and  the  physical  optics 
approximation . 
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I.  INTRODUCTION 


A.  HISTORY 

The  problem  of  scattering  from  periodic  rough  surfaces  has 
aroused  the  interest  of  physicists  and  engineers  for  the  last 
70  years  and  has  become  of  special  interest  in  recent  years, 
particularly  in  connection  with  the  propagation  of  radio  waves 
at  frequencies  above  30  MHz.  The  problem  of  scattering  from 
rough  surfaces  arises,  among  others,  in  line-of-sight  as  well 
as  mobile  radio  communication,  radar,  sea  scattering,  remote 
sensing,  tropospheric  propagation,  radio  astronomy,  and 
acoustics. 

Mathematically,  the  problem  of  the  scattering  of  waves  by 
a  rough  surface  is  difficult  to  formulate  and  even  more 
difficult  to  obtain  nvimerical  results. 

B.  ON-SURFACE  RADIATION  CONDITION  (OSRC)  METHOD 

Recently,  a  new  formulation  of  electromagnetic  scattering 

has  been  introduced  [Ref.  l].  This  formulation  of 

electromagnetic  wave  scattering  by  convex,  two-dimensional 
conducting  bodies,  the  On-Surface  Radiation  Condition  (OSRC) 
approach,  is  based  upon  an  expansion  of  the  radiation 
condition  applied  directly  on  the  surface  of  a  scatterer. 
According  to  the  OSRC  approach,  the  application  of  a  suitable 
radiation  condition  directly  on  the  conducting  scattering 
surface  can  lead  to  substantial  simplification  of  the 
frequency-domain  integral  equation  for  the  scattered  field. 
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For  the  transverse  magnetic  (TM)  case,  the  integrand  is  known 
explicitly.  For  the  transverse  electric  (TE)  case,  the 

integrand  can  be  constructed  by  solving  an  ordinary 
differential  equation  around  the  scatterer  surface  contour. 
Implicit  in  the  OSRC  formulation  is  the  fact  that  only 
outgoing  waves  are  asstimed  to  exist  on  the  surface  of  the 
scatterer.  Thus,  the  method  will  have  serious  drawbacks  in 
cases  where  the  latter  condition  is  not  met. 

C.  STATEMENT  OF  THE  PROBLEM 

In  this  thesis  we  apply  the  OSRC  technique  to  scattering 
by  a  periodic  rough  surface.  We  consider  a  plane 
electromagnetic  wave  illuminating  a  two-dimensional,  perfectly 
conducting  periodic  rough  surface  for  the  transverse  magnetic 
(TM)  polarization  case  (magnetic  field  parallel  to  the 
propagation  plane  or  electric  field  parallel  to  the  surface 
ridges  as  in  Figure  1)  . 

The  analysis  of  this  problem  includes  the  calculation  of 
the  induced  electric  current  as  well  as  the  echo  width  of  the 
scattering  surfaces.  Both  of  these  quantities  are  calculated 
in  Chapter  II. 

In  Chapter  III,  we  calculate  the  current  distribution  and 
the  echo  width  of  a  planar  thin  conducting  strip,  while  in 
Chapter  IV  we  calculate  the  same  two  quantities  for  a 
sinusoidal  scattering  surface  and  the  surface  current 
distribution  on  a  scatterer  of  triangular  shape. 
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The  numerical  results  of  the  current  distributions  as  well 
as  the  echo  width  of  the  planar  strip  and  the  sinusoidal 
surface  are  presented  in  Chapter  V.  Comparisons  of  these 
results  with  results  from  the  physical  optics  approximation 
and  with  the  exact  solution  (method  of  moments)  ,  are  also 
presented  in  this  chapter. 

In  Chapter  VI,  conclusions  and  recommendations  are 
presented.  Some  mathematical  derivations  and  the  computer 
codes  developed  to  evaluate  and  test  our  results  are  listed  in 
the  appendices. 
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IZ.  FORMULATION 


In  this  thesis  we  apply  the  OSRC  technique  to  scattering 
by  a  rough  surface.  In  this  chapter  we  will  formulate  the 
problem  and  calculate  both  the  surface  current  and  the  far- 
field  echo  width  of  a  periodic  rough,  perfectly  conducting 
surface. 

A.  NOMENCLATURE 

We  will  consider  a  plane  electromagnetic  wave  illuminating 
a  two-dimensional,  perfectly  conducting,  rough  surface  for  the 
transverse  magnetic  (TM)  polarization  case  (Figure  1) .  The 
width  of  the  rough  surface  is  considered  finite  and  edge 
effects  are  ignored. 

We  shall  use  Cartesian  coordinates  x,  y,  z  with  origin  0 
and  unit  vectors  St,  2.  The  two  dimensional  rough  surface 
is  specified  by  the  function  z  =  f(x)  where  the  mean  level  of 
this  surface  is  the  plane  z  =  0. 

All  quantities,  except  angles,  associated  with  the 
incident  field  will  be  denoted  by  the  superscript  i  and  those 
associated  with  the  sca,.tered  field  by  the  superscript  s. 
Thus,  the  incident  field  is  and  the  scattered  field  E®. 
The  medium  above  the  surface  is  assumed  to  be  free  space.  We 
shall  assume  to  be  linearly  polarized.  This  implies  that 

we  need  only  the  scalar  value  17®  of  the  vector  E^ . 
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B.  GENERAL  THEORY 

For  TM  polarization,  the  incident  electric  field 

propagating  perpendicular  to  the  x-z  plane  at  an  angle  0q 
with  respect  to  the  +z  axis,  as  shown  in  Figure  1,  is  given 
[Ref.  2] 

(2.1) 

where  is  the  magnitude  of  the  incident  electric  field 


-  zCX)^^) 

/C(j  is  the  wavenumber  of  the  free  space  defined  as 


(2.2) 


(2.3) 


and  harmonic  time  dependence  is  assumed.  The  quantity  Iq  is 
the  free  space  wavelength  associated  with  the  incident 
electromagnetic  wave  of  frequency  f^,  c  is  the  speed  of  light 
in  free  space,  and  o  =  is  the  radian  frequency  of  the 
incident  wave. 

The  scattered  electric  field  is  [Ref.  2] 


(2.4) 


where  U°  is  given,  in  the  far-field,  by  [Ref.  1] 


U”  = 


[SiCoXr] 


2  C 


dn' 


+  jiCgCOS6i/^ 


(2.5) 
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Figure  1  Scattering  geometry  for  a  TM  polarized  incident 
plane  wave. 


where  C  represents  the  contour  of  the  surface,  /i'  is  the  unit 
vector  normal  to  the  curve  C  at  point  P(x',z'),  as  shown  in 
Figure  2,  and  d  /dn'  denotes  an  outward  normal  derivative  on 
C.  Note  the  integration  in  Equation  (2.5)  is  with  respect  to 
the  arc  length,  s',  corresponding  to  the  point  P(x',z')  at 
distance  P  from  the  origin,  and  that  r  is  the  distance  of  the 
observation  point  from  the  origin.  The  unit  vector  f  in  the 
direction  of  the  observation  point  is 

f  =  sin0-)?  +  cos0-^  (2.6) 

where  0  is  the  scattering  angle,  measured  in  opposite  sense  to 
the  incidence  angle  0g. 
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Figure  2  Notation  for  various  vectors  and  angles. 


In  Equation  (2.5),  the  terms  cosS  and  Rcosf  are  given  by 

cos6  =  f-fi'  (2.7) 


i?cosi|r  =  R'f  (2.8) 

where 

R  =  x'Ji  +  z'2  (2.9) 
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The  unknown  quantity  dV^/dn'  can  be  expressed  in  terms  of 
the  incident  field  using  the  second  order  radiation  boundary 
condition.  It  is  given  by  [Ref.  l] 


3(7^  _[  C  ix')  ,  -r,  ,  l^ri^ _ 1 _ 

dn' ~  2  °  8  [-jJCo-C(xO  ]  J  2  [-jJCo-C(xO  ]  ds'^ 


(2.10) 


where  f(x')  is  the  curvature  of  the  two  dimensional 
scattering  surface  at  point  P(x',z')  and  3^  /ds*^  is  the 
second  derivative  with  respect  to  the  arc  length  of  C. 

The  curvature  f (x)  of  the  curve  C  that  is  defined  by  an 
equation  of  the  type  z  =  Jf(x)  is  given  by  [Ref.  3] 


while  the  quantity  /ds*^  is  shown  in  Appendix  A  to  be 


ds'^ 


-ic|sin20o  +  k^-^cos0oj2sineo  -  -^cos0oj  +  jk^-^cosB^ 


(sin0o  -  ^cos0o)-^4l 

_ °  dx  °  dx  dx^ 

L  .  I  dz\^]^ 


m] 


(2.12) 


By  substituting  Equations  (2.11),  (2.12),  into  Equation 

(2.10)  and  carrying  out  some  manipulation,  we  get 
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I 


du‘ 

dn' 


1  + 


d^z 

dx^ 

-I  -  _ 

d}z 

dx^ 

/ 

*  1 

\  1  dz\^' 

■§  d^z 

idx)  J 

4 

Jko 

K 

(dx)  . 

+  - 

dx^ , 

+  jK 


2 

2 


COS0, 


dz 

dx\ 


(2sin0f 


^cose.)) 


jJCc 


COS0„ 


1  + 


dz\^ 

dx) 


(2.13) 


The  quantity  dlf /dn*  is  expressed  in  terms  of  and  hence  is 
completely  )cnown  at  all  points  on  the  rough  surface. 

The  unit  vector  fi  * ,  normal  to  the  curve  C  at  point 
P(x',2'),  is  given  by  [Ref.  3] 


A  = 


/ 


± 


V 


(2.14) 


where  the  upper  sign  corresponds  to  the  outward  normal  and  the 
lower  sign  to  the  inward  normal. 

The  formula  for  the  partial  derivative  of  the  incident 
electric  field  with  respect  to  the  outward  normal  vector, 
which  we  will  use  later,  is  given  by 
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^  =  jku^ 

dn'  ^  ' 


(2.15) 


-^sinSg  +  COS0O 


Substituting  Equations  (2.6)  and  (2.14)  into  Equation 
(2.7),  we  get 


cos6  =  f‘fi  = 


--^sin0  +  COS0 
dx 


(2.16) 


Likewise  substituting  Equations  (2.9)  and  (2.6)  into 
Equation  (2.8),  we  get 

k^RcosT^  =  koiR’f)  =  iC(j(xsir20  +  zcosd)  (2.17) 

We  will  normalize  all  dimensions  to  the  free  space 
wavenumber  k^.  Accordingly,  we  let 

P  =  JCgX  ,  T  =  kgZ  (2.18) 

so  that 

dx  =-^dp  ,  dz  =  -^dz  (2.19) 

and 


dz 

dx 


dp 


cPz  _  .  d^T 
°dp2 


dx^ 


(2.20) 


Combining  Equations  (2.18),  (2.19)  and  (2.20)  with 
Equations  (2.2),  (2.16),  (2.17)  and  (2.13),  we  express  all 
quantities  in  terms  of  the  normalized  variables  p,T  as 
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Iji  =  -  »coaeo) 


— ^sin0  +  COS0 
cos5  =  — - 


kgRcos^  =  psin0  +  tcos0 


d^T 


dp' 


COS0„ 


The  arc  length  s  is  related  to  the  surface 


3] 


(2.21) 


(2.22) 


(2.23) 


+  J 


NIX 

(2.24) 

>  y 


z(x)  by  [Ref. 
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(2.25) 


SO  that  the  differential  arc  length  ds'  can  be  expressed  in 
terms  of  the  normalized  coordinates  as 


(2.26) 


Substituting  Equation  (2.26)  in  Equation  (2.5),  we  get: 


[SiCoiir] 


To  organize  our  work  better  and  to  enable  more  efficient 
structural  programming,  we  rewrite  the  above  equation  as 


U'  «  V(T)-f(F(p)-p(p))df> 


(2.28) 


Where 


V(r)  = 


-  i) 


(2.29) 


F(p)  =  +  jkgcosbu^ 

an' 


(2.30) 
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C)(p)  = 


jkoRcoB^ 


1  + 


(f) 


(2.31) 


and  duV^n^  ,  cosS,  u’  and  k^Rcos^  as  given  by  Equations 
(2.24),  (2.22),  (2.21)  and  (2.23)  respectively. 


C.  SURFACE  CURRENT  AND  ECHO  WIDTH 

Two  quantities  needed  to  compare  our  results  with  those  in 
the  literature  are  the  electric  surface  current  induced  on  the 
scattering  surface  and  the  far  field  echo  width  of  this 
surface.  These  quantities  will  be  presented  next. 

Since  =  U*  +  U*  and  the  total  magnetic  field  is  given 

by 


(2.32) 


the  normalized  electric  current  on  the  scattering  surface, 
equal  to  the  ratio  of  the  total  to  the  incident  magnetic 
field,  is 


—  =  1  ( du"  ^ 

^  dn  dn  ) 


(2.33) 


In  the  Equation  (2.33)  the  only  previously  undefined  term 
is  the  magnitude  of  the  incident  magnetic  field,  h’,  which  can 
be  obtained  in  terms  of  the  electric  field  using  the  Maxwell's 
equation  [Ref.  4] 


=  — ^(Vx^ 


wPoV 


dz  dx 


(2.34) 
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For  the  case  of  a  plane  incident  wave 


_  ■^0 


no 


(2.35) 


where  t7o  is  the  intrinsic  impedance  of  the  free  space. 

The  normalized  surface  electric  current  is  now  given  by 


1  (  dU' 


(2.36) 


The  far-field  echo  width  a  of  the  scattering  surface  is 
defined  as  [Ref.  2] 


o  =  lim  1 .1^.11  ■  2iir| 

1 1^1"  J 

Since  all  quantities  have  been  previously  defined, 
easily  calculated. 


(2.37) 

a  can  be 
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ZII.  THIN  CONDUCTING  STRIP 

The  inatheinatical  derivations  of  the  previous  chapter  are 
for  an  arbitrary  rough  surface.  In  this  chapter  we  specialize 
to  the  case  of  TM  scattering  by  a  thin  conducting  strip. 
Because  results  for  the  induced  current  and  echo  width  for  a 
conducting  strip  are  known  by  other  means  this  will  serve  as 
a  check  to  our  formulation. 

A.  SURFACE  CURRENT 

For  a  thin  planar  perfectly  conducting  strip  with  the 
geometry  of  Figure  3,  the  surface  is  described  by 

z  ~  ax  ,  X  =  ap  (3.1) 

where  a  is  the  slope  of  the  surface. 

Using  the  results  of  the  previous  chapter,  the  normalized 
surface  electric  current  density  on  this  surface  is  given  by 

^  asinVcoseoj 

Notice  that  the  current  distribution  is  )cnown  analytically. 
This  is  in  contrast  to  other  methods  where  differential  or 
integral  equations  must  be  solved  to  obtain  the  current 
distribution. 

Let  us  now  derive  a  formula  for  the  scattered  echo  width. 
According  to  Appendix  B,  we  can  express  the  integral  in 


^  L  sin^So  ♦acose(,(2sin6o-acos6o) 
*  [  2(1  ♦a*) 
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Flgur*  3  Planar,  thin  conducting  strip  profile, 
illuminated  by  a  TM  polarized,  plane  wave. 


Note  that  no  special  consideration  for  the  fields  is  given  at 
the  edges  x  «  i  L. 
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other  quantities  needed  to  evaluate  the  integral  are 


jji  _  g-iP(«in*o  -  ^cobOo) 


(3.5) 


k^Rcos'^  =  p(sin6  +  acos0) 


(3.6) 


1  dU‘ 
dn 


(3.7) 


where 


D  =  Bin^0Q  +  acos0(,(2sin0o  -  acosB^)  (3.8) 

Substituting  the  above  equations  into  Equation  (3.4),  we 

get 

C/*  »  V{r)2j\l - [1  +  a2]2  f  e^®Pdp  (3.9) 

I  2[1  ♦  a^]) 

where 


G  *  (sin0  -  sin0o)  +  a(cos0  +  cos0o)  (3.10) 

The  integral  in  Equation  (3.9)  can  be  evaluated  in  closed 
form.  Carrying  out  the  integration  and  substituting  for  V(r) , 
we  get 


U* 


4J 


1  ♦  a- 


SiCgiir 


sin(/c5LG) 

k„LG 
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B.  FAR-FIELD  ECHO  WIDTH 

Now  that  the  scattered  field  is  known,  the  echo  width  can 
be  easily  calculated  using  Equation  (2.37).  Substituting 
Equations  (2.2)  and  (3.11)  into  Equation  (2.37),  and 
normalizing  the  echo  width  a  to  the  free  space  wavelength  Ig, 
we  obtain 


-  2JCol|l  - 


koLG 


(3.12) 


The  above  formula  is  for  the  bistatic  echo  width  (receiver 
and  transmitter  are  considered  at  different  locations) . 

The  monostatic  case  (transmitter  and  receiver  at  the  same 
location)  is  a  special  case  of  the  bistatic  case,  where  the 
angle  of  incidence  is  related  with  the  observation  angle 
0  by 


0  =  2m  -  Go 


(3.13) 
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IV.  SCATTERING  FROM  TRIANGULAR  AND  SINUSOIDAL 
PERFECTLY  CONDUCTING  SURFACE 

In  this  chapter  we  investigate  TM  scattering  from 
perfectly  conducting  triangular  and  sinusoidal  surfaces. 
Expressions  for  the  electric  current  induced  for  both  the 
surfaces  and  the  echo  width  for  the  sinusoidal  case  alone  are 
developed.  The  sinusoidal  surface  is  generally  considered  as 
a  prototype  surface  for  studying  scattering  from  more 
complicated,  rough  surfaces.  The  scattering  surface  profiles 
are  shown  in  Figure  4. 

A.  TRIANGULAR  CASE 

For  the  case  of  the  triangular  profile,  as  shown  in  Figure 
4(a),  the  surface  is  described  by 

z{x)  =  h  -  ,  0  i  X  s  —  (4.1a) 

u  2 

z(x)  =  -3h  +  .  -|  <  X  i  cf  (4.1b) 

where  d  is  the  wavelength  of  the  periodical  surface.  The 
expression  for  the  normalized  surface  electric  current  density 
(Equation  (2.36))  is  simplified  to 
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Figure  4  Profiles  of  periodic  surfaces;  (a)  Triangular; 
(b)  Sinusoidal. 
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H 


Tot 


=  U^\ 


L\ 


sin^0o  +  -^coseo|2sin0o  - 

^cose,)' 

2 

1 .  (42  r 

\dxl  \ 

> 

where 


-^sin6o  +  cosOq 


(4.2) 


dz  _  _  4h 
dx  d 


0  i.  X  i.  ■$ 
2 


(4.3a) 


dz  _  _4h 
dx  d 


4  <  X  i  d 
2 


(4.3b) 


B.  SINUSOIDAL  CASE 

The  sinusoidal  scattering  surface  profile  (Fig.  4(b))  is 
described  by 

z(x)  =  hcos(wx)  (4.4) 


where  w  =  2ir/d  and  d  is  the  wavelength  of  the  sinusoidal 
surface.  Hence 


=  -hwsiniwx) 
dx 


- =  =  -hw*COS(W3c) 

dx^ 


(4.5) 


dx 

dp 


(4.6) 
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and 


k^z  =  k(^hcos{wx) 


(4.7) 


These  quantities  are  used  to  express  the  normalized 
surface  current  induced  on  the  surface  and  defined  in  Equation 
(2.36).  The  echo  width  of  this  surface  is  obtained  in  the 
following  fashion. 

Since  the  scattering  surface  is  perfectly  conducting  with 
considerable  thickness,  the  integral  in  Equation  (2.28)  is 
shown  in  Appendix  B  to  be 

pm*y-  r 

'0{p)dp  *  J 

c  p— 

The  far-zone,  scattered  field  reduces  to 


(4.8) 


U‘  = 


p-**oif 

V(r)  f 


P— iCoU  ^0 


dU* 

dn 


1  + 


m] 


where 


kfjRcosr\f  =  psin9  +  kg/icos^-^jcose 


(4.10) 


The  incident  field  on  the  surface  of  the  scatterer  is 


-j(pBin6o  ~  JcgJ>coa|.Sjco8eoj 


(4.11) 


For  computational  purposes  it  is  convenient  to  rewrite  the 
incident  field  as 
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(4.12) 


where 


r  =  p  (sin6  -  ain0o) 


(COS0  -  COS0o) 


(4.13) 


Substituting  these  in  Equation  (4.9),  we  obtain 

p-*ktL 

U'  ^  V(T)  j*  if(p)dp  (4.14) 

p— M- 

where 


(4.15) 


We  normalize  the  limits  of  integration  to  the  quantity  Ic^’L. 
After  some  manipulation,  Equation  (4.14)  reduces  to 


U‘  =  V(r) 


*k^ 


J  H(p)  dp  *  j  Hip)  dp 

[-koL  0 


-  £7*  =  V(r)  f[H{p)  *H{-P)]dp  (4.16) 


With  the  substitution  p  =  k^Lu  ,  we  get 

V  =  V(r)kt,LM 

where 


(4.17) 


1 

j[HikoLu)  +  Hi-koLu)]du  (4.18) 

0 

Since  a  closed  form  formula  has  not  been  found  for  the 
scattered  field,  the  integration  in  Equation  (4.18)  is 
performed  numerically. 
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Once  the  quantity  M  has  been  calculated,  the  echo  width  a 
normalized  to  the  wavelength  is  obtained  from 

®  (4.19) 

Xo  Bit 

When  normalized  with  respect  to  the  length  (2L)  of  the 
scattering  surface,  the  echo  width  is 
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V.  NUMERICAL  RESULTS 


In  this  chapter,  computed  results  for  the  induced 
electric  surface  current  and  the  echo  width  are  presented. 
In  order  to  check  our  theory  and  computer  codes,  the  results 
of  the  analytical  solution  (OSRC)  are  compared  against  the 
moment  method,  treated  here  as  exact,  and  the  physical  optics 
solution.  Details  of  these  comparisons  are  presented  in  the 
following  sections. 

In  all  models  where  numerical  results  are  presented,  the 
following  factors  are  considered: 

1.  The  incident  electromagnetic  wave  is  a  Transverse 
Magnetic  plane  wave  incident  in  the  x~z  plane. 

2.  However,  in  an  actual  periodic  surface  the  surface  is 
infinite  in  the  direction  of  periodicity,  the  width  of  the 
scattering  surface  considered  here  is  finite  and  edge  effects 
are  ignored. 

3.  The  induced  electric  current  is  plotted  against  the 
normalized  distance  (across  one  period  of  the  periodic 
surface) . 

4.  The  resulting  echo  width  a  is  normalized  to  the 
wavelength  for  the  planar  strip  and  to  the  surface  total 
length  2L  for  the  sinusoidal  surface. 

5.  Curves  presented  for  the  exact  and  physical  optics 
solution  are  taken  from  [Ref.  5]  and  [Ref.  6]. 
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A.  SURFACE  CURREMT  COMPARISONS 

For  a  planar  strip  (Figure  3)  of  length  2L,  the  magnitude 
of  the  surface  current  has  a  uniform  distribution  while  the 
phase  drops  linearly  across  the  strip,  as  shown  in  Figure  5, 
for  incidence  angles  0  =  0°  ,30°,  60°,  90°.  Notice  that  the 
current  distribution  does  not  exhibit  any  special  behavior  at 
the  edges. 

The  magnitude  of  the  surface  current  as  a  function  of  the 
incidence  angle  6q  for  the  OSRC  solution  as  compared  with  the 
physical  optics  solution  is  shown  in  Figure  6.  We  can  see 
that  the  agreement  is  remarkable  between  the  two  techniques 
for  angles  of  incidence  close  to  zero.  However,  the  results 
deviate  for  angles  of  incidence  exceeding  45°.  For  normal 
incidence  the  numerical  (OSRC)  results  agreed  with  the  exact 
solution  to  an  accuracy  better  than  0.01  percent. 

The  comparisons  for  the  sinusoidal  surface  between  the 
numerical  results  generated  by  the  present  approach  with  both 
the  exact  solution  and  the  physical  optics  approximation  are 
shown  in  Figures  7,  8  and  9.  Data  pertaining  to  the  latter 
two  techniques  are  taken  from  [Ref.  5]. 

The  surface  currents  obtained  by  all  three  methods  for 
normal  incidence  are  shown  in  Figures  7  and  8.  Agreement  of 
our  solution  with  the  other  two  is  best  for  h  =  0.25A.g  and 
d  =  3.9l(,  (Figure  8),  while  for  h  =  0.25Xg  and  d  =  1.9A.g 
(Figure  7)  the  numerical  results  deviate  considerably  from  the 
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other  two.  The  results  seem  to  be  closer  to  the  physical 
optics  approximation  in  both  cases.  The  maximvim  deviation 
occurs  at  the  peak  and  trough  of  the  surface  where  the 
curvature  is  maximum.  This  is  expected  of  the  OSRC  technique. 

In  general  both  outgoing  and  incoming  waves  exist  on  the 
surface  of  the  scatterer.  Millar  [Ref.  7]  has  shown  that  for 
2iTh/d  <  0.448,  the  fields  on  the  surface  are  predominantly  of 
outgoing  type.  In  the  OSRC  technique,  a  higher  order 
radiation  boundary  condition  is  imposed  directly  on  the 
surface  of  the  scatterer.  This  implies  that  only  outgoing 
waves  have  been  assumed  to  exist  on  the  surface  of  the 
scatterer.  In  the  case  of  Figure  8,  2irh/d  =  0.403,  which  is 
less  than  the  critical  value  of  0.448.  Hence  better  agreement 
is  seen  for  this  case.  However,  in  Figure  7,  2irh/d  =  0.826 
resulting  in  poorer  agreement. 

In  Figure  9  we  present  the  results  for  the  same  geometry 
as  in  Figure  8,  but  for  a  75°  angle  of  incidence.  The  OSRC 
solution  is  not  affected  by  the  occurrence  of  shadowing  as  in 
the  case  of  the  physical  optics.  However,  the  results  differ 
from  the  exact  results  especially  in  the  central  region 
(shadow  region) . 

The  current  distributions  on  periodic  surfaces  with 
sinusoidal  and  triangular  profiles  (see  Figure  4)  are  shown  in 
Figures  10  through  15  represent.  Each  surface  has  a  period 
d  =  0.2Xq,  and  a  maximum  amplitude  h  =  O.IIq  (2irh/d  =  tt)  . 


27 


The  angles  of  incidence  considered  are  0^  =  0°,  30°  and  60°. 
In  each  figure,  the  OSRC  solution  results  are  shown  along  with 
the  exact  solution  as  given  in  [Ref.  6].  It  is  seen  from 
Figure  10  that  the  OSRC  technique  predicts  an  amplitude  of 
approximately  6,  both  at  the  crest  and  the  trough  of  a 
sinusoidal  profile,  whereas  the  exact  solution  predicts  values 
of  5  and  0  at  the  crest  and  trough,  respectively.  The  rather 
large  deviation  of  the  results  at  the  trough  is  attributed  to 
the  presence  of  standing  waves,  which  the  OSRC  technique  fails 
to  take  into  account. 

For  oblique  incidences,  OSRC  technique  predicts  a  jiimp 
discontinuity  in  the  current  magnitude  at  the  center  for 
triangular  profiles  as  observed  in  Figures  13  through  15. 
This  is  due  to  the  discontinuity  of  the  surface  slopes  at 
these  locations.  The  exact  solution  exhibits  the  proper  edge 
condition  whereas  the  OSRC  fails  to  do  so.  A  similar  jump 
discontinuity  is  predicted  by  the  physical  optics  [Ref.  6]. 

B.  ECHO  WIDTH  COMPUTATIONS  •  COMPARISONS 

The  normalized  monostatic  echo  width  a/k^  of  a  plane  strip 
with  length  2L,  as  a  function  of  the  incidence  angle  Bq  is 
shown  in  Figure  16.  The  normalized  bistatic  echo  width  c/Xq, 
as  a  function  of  the  observation  angle  0,  for  k^L  =  5,  k^L  = 
10  and  angle  of  incidence  0°  and  60°  is  shown  in  Figure  17. 

Figure  16  shows  also  the  comparison  between  the  numerical, 
the  exact  [Ref.  1],  and  the  physical  optics  [Ref.  8]  results. 
We  see  that  we  have  a  very  good  agreement  (within  1  dB)  with 
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the  physical  optics  solution  for  angle  of  incidence  up  to  70°. 
Agreement  with  the  exact  solution  is  poor  beyond  30°  angle  of 
incidence.  The  deviation  from  the  exact  solution,  for  larger 
angles  of  incidence,  is  probably  due  to  edge  currents  that 
have  been  ignored  in  the  formulation. 

The  normalized  monostatic  echo  width  a/2L,  of  a  sinusoidal 
surface  of  length  2L  with  parameters  h  =  0.25lo  and  d  =  S.OIq 
is  shown  in  Figure  18  while  Figure  19  shows  the  normalized 
bistatic  echo  width  a/2L  of  a  sinusoidal  surface  with  the  same 
profile,  for  angle  of  incidence  0°  and  60°.  No  comparison  is 
made  for  the  echo  width  of  the  sinusoidal  surface. 


29 


Normalized  Distance  x/L 


Normalized  Distance  x/L 


Figur*  5  Normalized  surface  current  versus  distance  ^cross 

a  horizontal  strip,  L  =  0.21q.  Lines:  ( )  “  ©o  ®  P  ' 

{-+-)  -  00  *  30";  (  X-)  -  00  «  60  ;  (••o--)  -  0o  =  90°.  (a) 

Magnitude;  (b)  Phase. 


Magnitude 


Figur*  6  Normalized  surface  current  versus  angle  of 
incidence  across  a  horizontal  thin  strip.  Solid  line: 
OSRC  results;  broken  line:  physical  optics  approximation. 
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0  0.2  0.4  0.6  0.8  1 

Normalized  Distance  x/d 
(b) 

7igur«  7  Normalized  surface  current^  versus  distance  across 
period  of  sinusoidal  surface,  Bq  =  0  ,  d  =  h  «  0.251q. 

Solid  lines:  OSRC  results;  broken  lines:  (-X-)  -  exact, 

(.  — o-  .)  -  physical  optics. 
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0  0.2  0.4  0.6  0.8  1 

Normalized  Distance  x/d 
(b) 

Figure  8  Normalized  surface  current^  versus  distance  across 
period  of  sinusoidal  surface,  =  o”,  d  =  3.9A,j,,  h  =  0.25A.Q. 
Solid  lines:  OSRC  results;  broken  lines:  (-X-)  -  exact, 

(  -o-  •)  -  physical  optics. 
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0  0.2  0.4 

Normalized  Distance  x/d 
(b) 

Figure  9  Normalized  surface  current  >gersus  distance  across 
period  of  sinusoidal  surface.  Bp  *  75  ,  d  =  S.SXg,  h  = 
O.ZSIq.  Solid  lines;  OSRC  results;  broken  lines:  (-x-)  - 
exact,  (•  -o-  •)  -  physical  optics. 


Normalized  Distance  x/d 
(b) 

Figure  10  Normalized  surface  current  versv^s  distance 
across  period  of  sinusoidal  surface,  0q  =  o“,  d  = 
h  =  O.IIq.  Solid  lines:  OSRC  results;  broken  lines:  exact 
solution. 


Normalized  Distance  x/ d 
(b) 

Figur*  11  Normalized  surface  current  versu^  distance 
across  period  of  sinusoidal  surface,  6q  -  30°,  d  =  0.2>.q, 
h  =  O.iXq.  Solid  lines:  OSRC  results;  broken  lines:  exact 
solution. 


Magnitude 


Normalized  Distance  x/d 
(b) 

Figura  13  Normalized  surface  current  versus  distance 
across  period  of  triangular  surface,  Bj,  *  0®,  d  *=  0.2Xo» 
h  »  O.lXg.  Solid  lines:  OSRC  results;  broken  lines:  exact 
solution. 
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Normalized  Distance  x/d 
(b) 

Figure  15  Normalized  surface  current  versus  distance 
across  period  of  triangular  surface,  6q  «  6o  ,  d  «  0.21q, 
h  =  O.iXq.  Solid  lines;  OSRC  results;  broken  lines:  exact 
solution. 
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Figur*  16  Normalized  monostatic  echo  width  o/1q,  of 
horizontal  conducting  strip.  Solid  lines;  OSRC  results 

broken  lines:  (•x-)  -  exact  solution,  ( - )  -  physical 

optics.  (a)  kflL  =  5;  (b)  k^L  «  10. 
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Figure  17  Normalized  blstatic  echo  width  o/lg  of 
horizontal  conducting  strip  -  Numerical  analysis  (OSRC) 
results.  Solid  lines:  k^L  «  10;  broken  lines:  k(,L  =  5. 
(a)  00=0°;  (b)  00=60°. 
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Figure  18  Nomalized  nonostatic  echo  width  a/2L  of 
einusoidal  surface,  d  -  3.91.,  h  -  0.25lo>*  “  OSRC  results 
Solid  line:  kgL  >  10;  broken  line:  k^L  «  5. 
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Observation  Angle  (degrees) 
(a) 


(b) 


7igur«  19  Nornallzed  blstatic  acho  width  o/2L  of 
sinusoidal  surface,  d  >  3.9Aq,  h  >  0.25io;  •  OSRC  results 
Solid  lines:  lc,,L  ■  10;  broken  lines:  k^L  «  5.  (a)  6o  =  o®; 

(b)  -  60®. 
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VI 


CONCLUSIONS 


A.  SUMMlOtY 

Electromagnetic  scattering  from  rough  surfaces  is  hard  to 
formulate  and  even  harder  to  solve.  In  the  On-Surface 
Radiation  Boundary  Condition  a  higher  order  radiation 
condition  is  imposed  directly  on  the  surface  of  the  scatterer. 
The  frequency-domain  integral  equation  for  the  scattered  field 
is  reduced  to  just  a  line  integral.  This  leads  to  substantial 
simplification  of  the  problem.  Analytical  expressions  are 
obtained  for  the  surface  scattered  fields. 

In  this  thesis,  the  OSRC  method  has  been  used  to  formulate 
electromagnetic  scattering  from  a  periodic  rough,  perfectly 
conducting  surface.  Transverse  magnetic  (TM)  polarized,  plane 
incident  electromagnetic  wave  was  considered  throughout  the 
whole  formulation;  the  surface  was  considered  of  finite 
length,  and  edge  effects  were  ignored. 

Three  different  scattering  surface  profiles  have  been 
considered  to  test  our  formulation.  The  surface  currents  as 
well  as  the  echo  width  of  these  surfaces  were  calculated  by 
developing  computer  codes  and  the  results  were  compared  with 
results  from  the  exact  solution  and  the  physical  optics 
approximation. 

For  all  the  examined  surfaces,  the  OSRC  method  predicts 
results  somewhat  between  the  two  other  solutions.  For  the 
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case  of  a  thin  strip,  the  OSRC  results  for  the  echo  width  are 
within  0.1  dB  (except  at  nulls)  of  the  exact  solution  for 
angles  of  incidence  between  0°  and  30°  and  deteriorating  for 
larger  angles  of  incidence.  Actually  OSRC  results  are  more 
similar  to  the  results  given  by  physical  optics. 

For  the  cases  of  a  triangular  surface  and  a  sinusoidal 
surface  with  parameters  h  and  d  such  that  2nh/d  >  0.448,  the 
OSRC  results  are  not  very  accurate. 

For  the  case  of  a  sinusoidal  surface,  the  validity  of  the 
OSRC  results  was  shown  to  depend  on  surface  profile  parameters 
as  well  as  the  angle  of  incidence.  When  the  quantity  2nh/d 
<  0.448  and  the  angle  of  incidence  is  close  to  zero,  the  OSRC 
results  are  very  accurate,  while  for  larger  angles  of 
incidence  it  is  somewhat  closer  to  the  exact  solution  than  the 
physical  optics  one.  Overall,  the  results  of  the  OSRC 
technique  seem  to  mimic  those  of  the  physical  optics  rather 
than  the  exact  solution.  This  is  clearly  seen  in  all 
comparisons  of  current  distribution  and  echo  width. 

B.  RECOMMENDATIONS 

With  the  basic  theory  for  the  TM  case  tested,  we  recommend 
the  investigation  of  surface  currents  and  echo  width  of  other 
periodic  surfaces.  The  formulation  and  the  validation  of  the 
same  problem,  for  the  Transverse  Electric  (TE)  case  is  another 
recommendation.  Finally,  how  much  the  edge  and  the  finite 
surface  length  affect  the  results  needs  to  be  investigated. 
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APPENDIX  A 


MATHEMATIC  DERIVATIONS 

In  this  Appendix  we  will  derive  the  second  partial 
derivative  of  the  incident  field  with  respect  to  the  arc 
length  on  the  scattering  surface. 

Let  the  electric  field  of  the  incident  wave  be  denoted  by 
an  arbitrary  function  of  the  type  U  =  U(x,z) .  The  surface  cf 
the  scatterer  is  described  by  z  =  z (x) . 

dU  _  ^  ^ . . 

ds  dx  ds  dz  ds  \  dx  dz  dxj  ds 

where  s  is  the  arc  length  on  the  scattering  surface. 

Since 


ds 

dx 


we  have  then 


dU  ^  dU  dz 
dU  dx  dz  dx 


If  we  name  W  =  Su/ds,  then 


_  dW 
ds^  ds 


dW  ^  dW  dz 
dx  dz  dx 


(A. 2) 


(A. 3) 


(A. 4) 
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It  can  be  shown  that 


dx 


d^U  ^  dz  ^  dU  <Pz\  IdU  ^  dU  dz\  dz  (Pz 

dx^  *  dxdz  dx  dzdx^l  \3x  dz  dx/  dx  dx^  (A. 5) 


and 


'  S^U  ^  d^U  dz" 
^  ^  [dzdx  dz^  dx, 


For  the  case  of  a  plane  incident  wave 


(A.6) 


_  Q-jk^^ixaitA^-zCO^^) 


(A.7) 


we  have 


dU- 

dx 


=  ■  *«*>  =  -jk^U^ 


{A.8) 


du^ 

dz 


-jJc^e =  jk,U^ 


(A. 9) 


^  -  *•")  =  -klu^ 

3x2 


(A.  10) 


^  =  kle-^(^^  -  *'*)  =  -klu^ 
3z2 


(A. 11) 
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(A.12) 


=  Jr  If  e'^C***  ■  *•*) 

dxdz  ■  dzdx 


where 


Jc,  =  iCoSinOo  ,  ic^  =  JcoCosS, 


(A.13) 


Finally,  substituting  Equations  (A. 5)  through  (A.13)  in 
Equation  (A. 4),  we  get 


ds^ 


= 


-J:osin^0(j  +  Jc|-^cos0oj2sin6o  -  -^cos0o|  +  jk^-^^cosQ 


-(tr 


+  jK 


(sin0o  -  -^cos0(,) 


dz  cPz 
dx  dx^ 


(A. 14) 
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APPENDIX  B 


EVALUATION  OF  INTEGRALS 

In  this  Appendix  we  present  an  approach  for  evaluating  the 
line  integrals  with  respect  to  the  arc  length  that  are  encountered 
in  Chapters  3  and  4  for  the  case  of  a  planar  strip  and  sinusoidal 
surface. 

We  are  interested  in  evaluating  an  integral  of  the  form 
/g(  )dp  where  C  is  a  conformal  contour  enclosing  the  surface. 

We  first  consider  the  planar  strip.  The  closed  contour  C 
around  the  scattering  surface  is  considered  to  comprise  of  four 
contours  C^,  Cj,  Cj  and  as  shown  in  Figure  Bl. 

The  contribution  from  contours  Cj  and  approaches  zero  for 
bounded  integrands.  The  arc  length  s,  as  well  as  the  differential 
arc  length,  ds,  on  the  two  contours  C,  and  Cg,  with  respect  to  the 
point  A  (Figure  Bl)  is 

s  =  ^+  x,ds  =  dx  (B.l) 

on  C, ,  and 

s  =  L  +  -  xj  ,  ds  =  -dx  (B.2) 

on  Cj. 

For  computational  purposes  let  us  name  the  integral  in 
Equation  (2.5)  as 

I  =  jiCoC/^COssjdS  (B.3) 
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Figure  Bl.  Integration  contour  for  a  thin  strip. 


We  then  have 

where  i?,  and  ^2  are  the  unit  normal  vectors  outward  to  the  contours 
C,  and  Cj  respectively. 

Using  Equations  (B.l)  and  (B.2),  we  see  that  Equation  (B.4) 
becomes 

Cl  Cj  W  } 

Since  Cj  =  -C,  and  -  -^2,  we  have 

Cl  -Cl  yCi  -Cl  j 


(B.5) 


(B.6) 
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or 


I 


Cl  Cl 


(B.7) 


Hence,  for  the  case  of  a  thin  conducting  strip,  we  have 

*  jkgU^coafAds  =  ij^^dx  (B.8) 

'  Cl 

Let  us  now  examine  the  case  where  the  thickness  of  the 
scattering  surface  is  no  longer  small,  so  that  the  second  contour 
Cj  is  considered  to  lie  inside  the  conducting  surface  (Figure  B2) . 
In  this  case  we  have 


J  a 


Ci*c,  \Ci*Ca  ^ 


(B.9) 


The  second  term  in  the  above  Equation  (B.9)  vanishes  on  the  closed 
contour  C^+Cj  as  shown  previously.  Inside  a  perfect  conductor, 
the  total  magnetic  field  is  zero.  Hence  on  the  contour  Cj,  d\J^/dn 
=  -du'/dn.  Substituting  these,  we  see  that  Equation  (B.9) 
becomes 


(B.IO) 


Since  C,  =  -Cj  and  =  -^2  -  it  follows  that 

C '  Cl  ^  Cl  2  Cj ' 
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APPENDIX  C 


PROGRAMS 

In  this  Appendix  we  present  the  computer  programs,  all  in 
FORTRAN  language  that  are  needed  to  perform  the  various 
calculations. 

A.  DISCUSSION 

Programs  STRIPCUR,  SINUSCUR  and  TRIANCUR  calculate  the 
surface  currents  induced  on  the  scattering  surfaces,  for  a 
strip,  sinusoidal  and  triangular  profiles,  respectively. 

Programs  STRIPMON  and  STRIPBI  calculate  the  normalized 
monostatic  and  bistatic  echo  width,  respectively,  of  thin, 
planar,  perfectly  conducting  strip.  Programs  SINUSMON  and 
SINUSBI  calculate  the  same  quantities  but  for  a  sinusoidal 
surface. 

Since  we  deal  with  complex  niimbers  most  of  the  above 
mentioned  programs  contain  subroutines  to  calculate  the  real 
(i.e.  STRPREAL,  SINREAL,  FRSTREAL  etc.)  and  the  imaginary  part 
(i.e.  STRPIMAG,  SINIMAG,  MOODIMAG  etc.)  of  a  quantity  and  use 
subroutine  PHASER  to  calculate  the  phase  of  this  complex 
quantity. 

Program  TRIANCUR  calculates  the  current  distribution  on 
the  triangular  surface,  by  examining  separably  the  two  halves 
of  the  surface,  the  first  with  negative  slope  (subroutines 
FRSTREAL,  FRSTIMAG) ,  and  the  second  with  positive  slope 
(subroutines  SCNDREAL,  SCNDIMAG) . 

Programs  SINUSMON  and  SINUSBI  carry  out  the  integration  of 
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Equation  (4.18),  using  Filon's  method,  as  it  is  given  by 
function  subprogram  FILON  [Ref.  9].  To  take  advantage  of  the 
fast  and  acurate  results  of  function  FILON,  we  had  to  separate 
the  even  and  the  odd  part  of  the  examined  function,  as  they 


are  shown 

in 

subroutines 

MOEVREAL, 

MOEVIMAG ,  BIEVREAL , 

BIEVIMAG, 

and 

MOODREAL, 

MOODIMAG, 

BIODREAL,  BIODIMAG 

respectively. 
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B.  PROGRAM  8TRZPCTJR 


THIS  PROGRAM  CALCULATES  THE  NORMALIZED  SURFACE  CURRENT 
ON  A  PLANAR,  PERFECTLY  CONDUCTING  STRIP,  DESCRIBED  BY 
THE  EQUATION:  2(x)  =  A  *  x 


INPUT  DATA: 

F  =  FREQUENCY  IN  GHz 

INCID  =  ELECTRIC  FIELD  ANGLE  OF  INCIDENCE  (IN  DEGREES) , 
WITH  RESPECT  TO  THE  VERTICAL  (+Z)  AXIS 

A  =  SLOPE  OF  THE  STRIP 

L  =  CONSIDERED  LENGTH  OF  THE  SURFACE 

M  =  NUMBER  OF  EQUIDISTANT  POINTS  BETWEEN  X=0  AND 

X=L  WHERE  OUR  FUNCTION  IS  GOING  TO  BE 
EVALUATED. 


OUTPUT  DATA: 

FREAL,FIMAG  *  REAL  AND  IMAGINARY  PART  RESPECTIVELY  OF 

THE  NORMALIZED  SURFACE  CURRENT 

AMAGNIT,  PHASE  =  MAGNITUDE  AND  PHASE  RESPECTIVELY  OF 

THE  NORMALIZED  SURFACE  CURRENT 


REAL  K,  INCID,  INCR,  L 

COMPLEX  ZET 

COMMON  K,  A,  THETAO 

OPEN  (UNIT=11,  FILE='STRIPCUR.DAT' ,  FORM= ' FORMATTED ' , 
&  STATUS* • UNKNOWN • ) 

READ(5,*)  F,INCID,A,L,M 

PI  =  4.*ATAN(1.) 

RAD  *  PI/180. 

THETAO  =  INCID*RAD 
K  =  20.*PI*F/3. 

INCR  =  L/M 

DO  100  X  =  0.,  L,  INCR 
DISTANCE  *  X/L 
FREAL  =  STRPREAL(X) 

FIMAG  =  STRPIMAG(X) 

ZET  =  CMPLX (FREAL, FIMAG) 

AMAGNIT  =  CABS (ZET) 
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PHASE  =  PHASER (FREAL,  FIMAG) 

WRITE (11, 20)  DISTANCE,  FREAL,  FIMAG,  AMAGNIT,  PHASE 
20  F0RMAT(1X,  F12.6,  4X,  P12.6,  4X,  F12.6,  4X,  F12.6, 

&  4X,  F12.6) 

100  CONTINUE 
STOP 
END 
C 

C  ****  It  It  *******11  hit -kit****  it  hit 

REAL  FUNCTION  STRPREAL(X) 

C  ************************* 

c 

REAL  K,  KC3 

COMPLEX  J,  POWER,  FUSCAT,  FUINC,  UI,  FUNC 
COMMON  K,  A,  THETAO 
C 

STHO  *  SIN (THETAO) 

CTHO  *  COS (THETAO) 

Z  *  A*X 

J  =  CMPLX(0.,1.) 

POWER  =  -  J*K*(X*STH0  -  Z*CTH0) 

C  *  SQRT(1.+  A*A) 

D  =  (STH0)**2.+  A*CTH0*(2.*STH0  -  A*CTH0) 

FUSCAT  =  J*(l.  -  D/(2.*C*C)) 

FUINC  =  J*((A*STH0  +  CTH0)/C) 

UI  *  CEXP( POWER) 

FUNC  =  -  J*UI* (FUSCAT  +  FUINC) 

FUNCREAL  =  REAL (FUNC) 

RETURN 

END 

C 

C  ************************* 

REAL  FUNCTION  STRPIMAG(X) 

C  ************************it 

c 

REAL  K,  KC3 

COMPLEX  J,  POWER,  FUSCAT,  FUINC,  UI,  FUNC 
COMMON  K,  A,  THETAO 
C 

STHO  =  SIN (THETAO) 

CTHO  *  COS (THETAO) 

Z  =  A*X 

J  =  CMPLX(0. ,1.) 

POWER  =  -  J*K*(X*STH0  -  Z*CTH0) 

C  =  SQRT(1.+  A*A) 

D  =  (STHO)* *2. +  A*CTH0*(2.*STH0  -  A*CTH0) 

FUSCAT  =  J*(l.  -  D/(2.*C*C)) 

FUINC  =  J*((A*STH0  +  CTH0)/C) 

UI  =  CEXP( POWER) 

FUNC  =  -  J*UI* (FUSCAT  +  FUINC) 

FUNCIMAG  =  AIMAG(FUNC) 

RETURN 
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END 


C  ************************** 

REAL  FUNCTION  PHASER (X,Z) 
************************** 


C 


c 


PI  «  4.*ATAN(1.) 


PI  =  ir 


IF  (X.EQ.O.O.AND.Z.EQ.0.0)  THEN 

ANGLE  =  0. 

ELSEIF  (X.EQ.O.O.AND.Z.GT.0.0)  THFN 

ANGLE  =  PI/2. 

ELSEIF  (X.EQ.O.O.AND.Z.LT.0.0)  THEN 

ANGLE  =  “PI/ 2. 

ELSEIF  (X.GT.O.O.AND.Z.GE.0.0)  THEN 

ANGLE  =  ATAN(ABS(Z/X) ) 
ELSEIF  (X.GT.O.O.AND.Z.LT.0.0)  THEN 

angle  =  -  ATAN(ABS(Z/X) ) 

ELSEIF  (X.LT.O.O.AND.Z.LT.0.0)  THEN 

angle  =  -PI/2.  -  ATAN(ABS(Z/X) ) 
ELSEIF  (X.LT.O.O.AND.Z.GE.0.0)  THEN 

ANGLE  =  PI  -  ATAN(ABS(Z/X) ) 

ENDIF 

COST  =  180/jr 

CONST  «  57.29577951 
PHASER  =  ANGLE* CONST 
RETURN 
END 
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C.  PROGRAM  SIMUSCUR 


THIS  PROGRAM  CALCULATES  THE  NORMALIZED  SURFACE  CURRENT 
ON  A  SINUSOIDAL,  PERFECTLY  CONDUCTING  SURFACE,  DESCRIBED 
BY  THE  EQUATION:  Z(X)  =  H*COS ( [2*pi/d] *X) 


INPUT  DATA: 

F  =  FREQUENCY  IN  GHz 

INCID  *  ELECTRIC  FIELD  ANGLE  OF  INCIDENCE  (IN  DEGREES) 
WITH  RESPECT  TO  THE  VERTICAL  (+Z)  AXIS 

H,  d  =  PARAMETERS  OF  THE  SINUSOIDAL  SURFACE 

M  =  NUMBER  OF  EQUIDISTANT  POINTS  BETWEEN  X=0  AND 

X=d  WHERE  OUR  FUNCTION  IS  GOING  TO  BE 
EVALUATED 


OUTPUT  DATA: 

FREAL,FIMAG  =  REAL  AND  IMAGINARY  PART  RESPECTIVELY  OF 
THE  NORMALIZED  SURFACE  CURRENT 

AMAGNIT,  PHASE  *  MAGNITUDE  AND  PHASE  RESPECTIVELY  OF 

THE  NORMALIZED  SURFACE  CURRENT 


REAL  K,  INCID,  INCR 

COPLEX  ZET 

COMMON  K,  H,  W,  THETAO 
C 

OPEN  (UNIT=12,  FILE=*SINUSCUR.DAT' ,  FORM= ' FORMATTED • , 
&  STATUS= • UNKNOWN ' ) 

READ(5,*)  F,  INCID,  H,  d,  M 
C 

PI  =  4.*ATAN(1.) 

W  =  2 . *PI/d 

RAD  =  PI/180. 

THETAO  =  INCID*RAD 

K  =  20.*PI*F/3. 

INCR  =  d/M 

DO  100  X  =  0. ,  d,  INCR 
DISTANCE  *  X/d 
FREAL  =  SINREAL(X) 

FIMAG  =  SINIMAG(X) 

ZET  =  CMPLX (FREAL,  FIMAG) 

AMAGNIT  =  CABS (ZET) 

PHASE  =  PHASER (FREAL,  FIMAG) 
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WRITE (12, 20)  DISTANCE,  FREAL,  FIMAG,  AMAGNIT,  PHASE 
20  F0RMAT(1X,  F12.6,  4X,  F12.6,  4X,  F12.6,  4X,  F12.6, 

&  4X,  F12.6) 

100  CONTINUE 
STOP 
END 
C 

REAL  FUNCTION  SINREAL(X) 

C  ************************ 

C 

REAL  K,  KC3 

COMPLEX  J,  POWER,  BJKC3,  DJE,  FUSCAT,  FUINC,  UI,  FUNC 
COMMON  K,  H,  W,  THETAO 
C 

STHO  =  SIN (THETAO) 

CTHO  *  COS (THETAO) 

WX  =  W*X 
SW  =  SIN(WX) 

CW  *  COS(WX) 

Z  *  H*CW 
J  *  CMPLX(0.,1.) 

POWER  =  -  J*K*(X*STH0  -  Z*CTH0) 

A  «  -H*W*SW 
B  «  -H*W*W*CW 
BA  «  ABS(B) 

C  «  SQRT(1.+  A*A) 

D  *  K*((STH0)**2.+  A*CTH0*(2.*STH0  -  A*CTH0) ) 

E  =  B*(CTH0  +  A* (STHO  -  A*CTH0) / (C*C) ) 

KC3  =  K*C*C*C 
BJKC3  »  BA  +  J*KC3 
DJE  =  D  +  J*E 

FUSCAT  =  BA/(2.*KC3)*(1  -  BA/ (4 . *BJKC3) )  +  J 
&  +  C/(2.*BJKC3)*DJE 

FUINC  =  J*((A*STH0  +  CTH0)/C) 

UI  =  CEXP( POWER) 

FUNC  =  -  J*UI* (FUSCAT  +  FUINC) 

FUNCREAL  *  REAL (FUNC) 

RETURN 

END 

C 

Q  ************************ 

REAL  FUNCTION  SINIMAG(X) 

C  ************************ 

c 

REAL  K,  KC3 

COMPLEX  J,  POWER,  BJKC3,  DJE,  FUSCAT,  FUINC,  UI,  FUNC 
COMMON  K,  H,  W,  THETAO 
C 

STHO  *  SIN (THETAO) 

CTHO  »  COS (THETAO) 

WX  =  W*X 
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SW  =  SIN(WX) 

CW  =  COS(WX) 

Z  =  H*CW 
J  =  CMPLX(0. ,1.) 

POWER  =  -  J*K*(X*STHO  -  Z*CTHO) 

A  =  -H*W*SW 
B  =  -H*W*W*CW 
BA  =  ABS(B) 

C  =  SQRT(1.+  A*A) 

D  =  K*( (STH0)**2.+  A*CTH0*(2.*STH0  -  A*CTHO) ) 
E  =  B*(CTHO  +  A*(STHO  -  A*CTHO)/ (C*C) ) 

KC3  =  K*C*C*C 
BJKC3  =  BA  +  J*KC3 
DJE  =  D  +  J*E 

FUSCAT  =  BA/(2,*KC3)*(1  -  BA/ (4. *BJKC3) )  +  J 
&  +  C/(2.*BJKC3)*DJE 

FUINC  =  J*((A*STHO  +  CTHO)/C) 

UI  =  CEXP( POWER) 

FUNC  =  -  J*UI* (FUSCAT  +  FUINC) 

FUNCIMAG  =  AIMAG(FUNC) 

RETURN 

END 
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PROGRAM  TRIANCUR 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


& 


c 


THIS  PROGRAM  CALCULATES  THE  NORMALIZED  SURFACE  CURRENT 
ON  A  TRIANGULAR,  PERFECTLY  CONDUCTING  SURFACE,  DESCRIBED 
BY  THE  EQUATION: 

z(x)  =  H  -  (4*H/d)*x  ,  for  0  <  x  <  d/2  (FIRST  HALF) 
z(x)  =  -3*H  +  (4*H/d)*x  ,  for  d/2  <  x  <  d  (SECOND  HALF) 


INPUT  DATA: 

F  =  FREQUENCY  IN  GHz 

INCID  =  ELECTRIC  FIELD  ANGLE  OF  INCIDENCE  (DEGREES) , 
WITH  RESPECT  TO  THE  VERTICAL  (+Z)  AXIS 

H,  d  =  PARAMETERS  OF  THE  TRIANGULAR  SURFACE 

M  =  NUMBER  OF  EQUIDISTANT  POINTS  BETWEEN  X=0  AND 

X=d  WHERE  OUR  FUNCTION  IS  GOING  TO  BE 
EVALUATED 


OUTPUT  DATA: 

FREAL,FIMAG  «  REAL  AND  IMAGINARY  PART  RESPECTIVELY  OF 

THE  NORMALIZED  SURFACE  CURRENT 

AMAGNIT,  PHASE  =  MAGNITUDE  AND  PHASE  RESPECTIVELY  OF 

THE  NORMALIZED  SURFACE  CURRENT 


REAL  K,  INCID,  INCR 

COMPLEX  ZET 

COMMON  K,  H,  d,  THETAO 

OPEN  (UNIT=14,  FILE= 'TRIANCUR. DAT ' ,  FORM= ' FORMATTED ' , 

STATUS= ' UNKNOWN ’ ) 

READ(5,*)  F,  INCID,  H,  d,  M 

PI  =  4.*ATAN(1.) 

RAD  =  PI/180. 

THETAO  =  INCID*RAD 
K  =  20.*PI*F/3. 

INCR  =  d/M 

DO  100  X  =  0. ,  d,  INCR 
DISTANCE  =  X/d 

IF  ( DISTANCE. GE. 0.0. AND. DISTANCE. LE. 0.5)  THEN 
FREAL  =  FRSTREAL(X) 

FIMAG  =  FRSTIMAG(X) 

ELSE 
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FREAL  =  SCNDREAL(X) 

FIMAG  =  SCNDIMAG(X) 

ENDIF 

ZET  =  CMPLX( FREAL, FIMAG) 

AMAGNIT  =  CABS (ZET) 

PHASE  =  PHASER (FREAL,  FIMAG) 

WRITE(14,20)  DISTANCE,  FREAL,  FIMAG,  AMAGNIT,  PHASE 
20  F0RMAT(1X,  F12.6,  4X,  F12.6,  4X,  F12.6,  4X,  F12.6, 

&  4X,  F12.6) 

100  CONTINUE 
STOP 
END 

■khliltlfkltlfkltlfkltltlflflfkltliltltlfklt 

REAL  FUNCTION  FRSTREAL(X) 

REAL  K,  KC3 

COMPLEX  J,  POWER,  FUSCAT,  FUINC,  UI,  FUNC 
COMMON  K,  H,  d,  THETAO 

STHO  =  SIN (THETAO) 

CTHO  =  COS (THETAO) 

A  =  -  4.*H/d 
C  =  SQRT(1.+  A*A) 

Z  =  H  +  A*X 
J  =  CMPLX(0.,1.) 

DI  *  (STH0)**2.+  A*CTH0*(2.*STH0  -  A*CTH0) 

FUSCAT  =  J*(l.  -  DI/(2.*C*C)) 

FUINC  =  J*((A*STH0  +  CTH0)/C) 

POWER  =  -  J*K*(X*STH0  -  Z*CTH0) 

UI  =  CEXP(IOWER) 

FUNC  =  -  J*UI* (FUSCAT  +  FUINC) 

FRSTREAL  =  REAL (FUNC) 

RETURN 
END 

************************* 

REAL  FUNCTION  FRSTIMAG(X) 

************************* 

REAL  K,  KC3 

COMPLEX  J,  POWER,  FUSCAT,  FUINC,  UI ,  FUNC 
COMMON  K,  H,  d,  THETAO 

STHO  =  SIN (THETAO) 

CTHO  =  COS (THETAO) 

A  =  -  4.*H/d 
C  =  SQRT(1.+  A*A) 

Z  =  H  +  A*X 
J  =  CMPLX(0. , 1. ) 

DI  =  (STH0)**2.+  A*CTH0*(2.*STH0  -  A*CTH0) 
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FUSCAT  =  J*(l.  -  DI/(2.*C*C)) 

FUINC  =  J*((A*STHO  +  CTHO)/C) 

POWER  =  -  J*K*(X*STHO  -  Z*CTHO) 

UI  =  CEXP( POWER) 

FUNC  =  -  J*UI* (FUSCAT  +  FUINC) 

FRSTIMAG  =  AIMAG(FUNC) 

RETURN 

END 

************************* 

REAL  FUNCTION  SCNDREAL(X) 

************************* 

REAL  K,  KC3 

COMPLEX  J,  POWER,  FUSCAT,  FUINC,  UI,  FUNC 
COMMON  K,  H,  d,  THETAO 

STHO  =  SIN (THETAO) 

CTHO  =  COS (THETAO) 

A  =  +  4.*H/d 
C  =  SQRT(1.+  A*A) 

Z  =  -3*H  +  A*X 
J  =  CMPLX(0. ,1. ) 

DI  *=  (STH0)**2.+  A*CTH0*(2.*STH0  -  A*CTHO) 
FUSCAT  =  J*(l.  -  DI/(2.*C*C)) 

FUINC  =  J*((A*STHO  +  CTHO)/C) 

POWER  =  -  J*K*(X*STHO  -  Z*CTHO) 

UI  =  CEXP( POWER) 

FUNC  *  -  J*UI* (FUSCAT  +  FUINC) 

SCNDREAL  =  REAL (FUNC) 

RETURN 

END 

************************* 

REAL  FUNCTION  SCNDIMAG(X) 

************************* 

REAL  K,  KC3 

COMPLEX  J,  POWER,  FUSCAT,  FUINC,  UI,  FUNC 
COMMON  K,  H,  d,  THETAO 

STHO  =  SIN (THETAO) 

CTHO  =  COS (THETAO) 

A  =  +  4.*H/d 
C  =  SQRT(1.+  A*A) 

Z  =  “3*H  +  A*X 
J  =  CMPLX(0. ,1.) 

DI  =  (STH0)**2.+  A*CTH0*(2.*STH0  -  A*CTHO) 
FUSCAT  =  J*(l.  -  DI/(2.*C*C)) 

FUINC  =  J*((A*STHO  +  CTHO)/C) 

POWER  =  -  J*K*(X*STHO  -  Z*CTHO) 

UI  =  CEXP( POWER) 
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FUNC  =  -  J*UI*(FUSCAT  +  FUINC) 
SCNDIMAG  =  AIMAG(FUNC) 

RETURN 

END 
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E.  PROGRAM  8TRIPM0M 


C  THIS  PROGRAM  CALCULATES  THE  MONOSTATIC  ECHO  WIDTH 

C  OF  A  THIN,  PLANAL,  CONDUCTING  STRIP  OF  LENGTH  2*L. 

C 
C 

C  INPUT  DATA: 

C 

C  F 

C  KOL 

C 

C  INCLIN 

C 

C  STEP 

C 
C 

C  OUTPUT  DATA: 

C 

C  ECHOMONO  *  NORMALIZED  -  WITH  RESPECT  TO  THE  WAVELENGTH  - 
C  MONOSTATIC  ECHO  WIDTH  AS  A  FUNCTION  OF  THE 

C  INCIDENT  ANGLE  THETAO,  WHICH  HERE  (MONOSTATIC 

C  CASE)  IS  EQUAL  TO:  2*PI  -  ANGLE  OF  INCIDENCE. 

C 
C 

REAL  KOL, INCLIN 

OPEN  (UNIT*1,  FILE=’STRIPMON.DAT' ,  FORM* ' FORMATTED ' , 

&  STATUS* • UNKNOWN ’ ) 

READ(5,*)  F, KOL, INCLIN, STEP 

C 

PI  =  4 . *ATAN ( 1 . ) 

RAD  =  PI/180. 

A  =  TAN (INCLIN) 

C2  =  1.  +  A*A 
C 

DO  100  QI  *  0.,  90.,  STEP 
THETAO  =  QI*RAD 
THETAR  =  2.* PI  -  THETAO 
STHO  =  SIN (THETAO) 

CTHO  =  COS (THETAO) 

STH  *  SIN (THETAR) 

CTH  =  COS (THETAR) 

C 

D  =  (STH0)**2.  +  A*CTH0*(2.*STH0  -  A*CTH0) 

G  =  STH  -  STHO  +  ALPHA* (CTH  +  CTHO) 

C 

ECHOMO=C2/ (2 . *PI) * (2 . * (1 . -D/ (2 . *C2 ) ) *SIN (K0L*G) /G) **2 
WRITE (1,*)  QI,ECHOMO 
100  CONTINUE 
STOP 
END 


=  FREQUENCY  IN  GHz 

=  K0*L  IS  THE  PRODUCT  OF  THE  STRIP  HALF 
LENGTH  (L)  TIMES  KO  (K0=2*PI/WAVELENGTH) 

=  STRIP  INCLINATION  (IN  DEGREES)  WITH  RESPECT 
TO  THE  HORIZONTAL  (+X)  AXIS 
=  ANGLE  OF  INCIDENSE  INCREMENT 
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PROGRAM  8TRIPBI 


THIS  PROGRAM  CALCULATES  THE  BISTATIC  ECHO  WIDTH 
OF  A  THIN,  PLANAL,  CONDUCTING  STRIP  OF  LENGTH  2*L. 


INPUT  DATA: 

F  =  FREQUENCY  IN  GHz 

KOL  =  KO*L  IS  THE  PRODUCT  OF  THE  STRIP  HALF 

LENGTH  (L)  TIMES  KO  (K0=2*PI/WAVELENGTH) 
INCLIN  =  STRIP  INCLINATION  (IN  DEGREES)  WITH  RESPECT 
TO  THE  HORIZONTAL  (+X)  AXIS 

INCID  =  ELECTRIC  FIELD  ANGLE  OF  INCIDENCE  (IN  DEGREES) 
WITH  RESPECT  TO  THE  VERTICAL  (+Z)  AXIS 
STEP  =  OBSERVATION  ANGLE  INCREMENT 


OUTPUT  DATA; 

ECHOBI  =  NORMALIZED  BISTATIC  ECHO  WIDTH  AS  A 

FUNCTION  OF  THE  OBSERVATION  ANGLE,  THETA, 
FOR  A  GIVEN  INCIDENT  ANGLE 


REAL  KOL, INCLIN , INCID 

OPEN  (tJNIT=2 ,  FILE=  *  STRIPBI .  DAT ' ,  FORM= '  FORMATTED '  , 

&  STATUS= ' UNKNOWN • ) 

READ(5,*)  F, KOL, INCLIN, INCID, STEP 

C 

PI  =  4.*ATAN(1.) 

RAD  =  PI/180. 

THETAO  =  INCID*RAD 
A  =  TAN (INCLIN) 

C2  =  1.  +  A*A 

C 

DO  100  QI  =  -180,  180,  STEP 
THETAD  =  QI 
THETAR  =  QI*RAD 
STHO  =  SIN (THETAO) 

CTHO  =  COS (THETAO) 

STH  =  SIN ^THETAR) 

CTH  =  COS (THETAR) 

C 

D  =  (STH0)**2.  +  A*CTH0*(2.*STH0  -  A*CTH0) 

G  =  STH  -  STHO  +  ALPHA* (CTH  +  CTHO) 

C 

ECHOBI=C2/ ( 2 . *PI ) * ( 2 . * ( 1 . -D/ ( 2 . *C2 ) ) *SIN (K0L*G) /G) **  2 
WRITE (2,*)  THETAD, ECHOBI 
100  CONTINUE 
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STOP 

END 
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6.  PROGRAM  SIMUSMOM 


THIS  PROGRAM  CALCULATES  THE  MONOSTATIC  ECHO  WIDTH  OF  A 
SINUSOIDAL  SURFACE  z(x)  =  H*COS ( [2*pi/d] *x) 


INPUT  DATA: 

FREQ  =  FREQUENCY  IN  GHz 

H,  d  =  PARAMETERS  OF  THE  SINUSOIDAL  SURFACE 

KOL  =  LIMIT  OF  INTEGRATION  (KOL  =  K*L) ,  WHERE  : 
"K”  IS  THE  PROPAGATION  CONSTANT  AND 
"L"  IS  THE  ON  THE  X-AXIS  SURFACE  PROJECTION 
LENGTH 

STEP  =  ANGLE  OF  INCIDENSE  INCREMENT 

M  =  NUMBER  OF  EQUIDISTANT  POINTS  BETWEEN  THE  TWO 
NORMALIZED  LIMITS  (0  AND  1)  OF  INTEGRATION 
(M  ODD,  M  <  999) . 


OUTPUT  DATA: 

ECHOMO  *  NORMALIZED  (WITH  RESPECT  TO  THE  WAVELENGTH) 
MONOSTATIC  ECHO  WIDTH,  AS  A  FUNCTION  OF  THE 
INCIDENT  ANGLE  (THETAO)  WHICH  HERE  (MONOSTATIC 
CASE)  IS  EQUAL  TO:  2*PI  -  OBSERVATION  ANGLE. 

ECHO  =  NORMALIZED  (WITH  RESPECT  TO  THE  SURFACE  LENGTH 
2*L)  MONOSTATIC  ECHO  WIDTH  AS  A  FUNCTION  OF  THE 
ANGLE  OF  INCIDENCE,  THETAO. 


REAL  F(999) 

REAL  K,K0L,K0LU,INCR,INCID,MAGN2_R 
COMPLEX  J,  R 

EXTERNAL  FILON,  FREAL,  FIMAG 
COMMON  K,  H,  W,  THETAO,  THETAR 

OPEN  (UNIT=7,  FILE='SINUSMON.DAT' ,  FORM= ' FORMATTED ' , 
&  STATUS= ' UNKNOWN ’ ) 

READ(5,*)  FREQ,  H,  d,  KOL,  STEP.  M 
C 

PI  =  4.*ATAN(1.) 

W  =  2.*PI/d 
RAD  =  PI/180. 

K  =  20.*PI*FREQ/3. 

INCR  =  l./M 


69 


o  o 


c 


DO  100  QI  =  0.,  90.,  STEP 
THETAO  =  QI*RAD 
THETAR  =  2.*PI  -  THETAO 
TSO  «  KOL*  (SIN  (THETAR)  -•  SIN  (THETAO)) 

SGN  =  SIGN (1, TSO) 

ABSTSO  »  ABS(TSO) 

C 

1  =  1 

DO  40  U  =  0. ,  1. ,  INCR 
KOLU  =  K0L*U 
F(I)  =  MOEVREAL(KOUJ) 

1  =  1  +  1 
40  CONTINUE 

R1  =  FILON(F, ABSTSO, 0. ,1. ,M, 999,1) 

C 

1  =  1 

DO  45  U  =  0. ,  1. ,  INCR 
KOLU  =  K0L*U 
F(I)  =  MOEVIMAG(KOLU) 

1  =  1  +  1 

45  CONTINUE 

R2  =  FILON (F, ABSTSO, 0. ,1. ,M, 999,1) 

C 

1  =  1 

DO  50  U  =  0.,  1.,  INCR 
KOLU  =  K0L*U 
F(I)  =  MOODREAL(KOLU) 

1  =  1  +  1 

50  CONTINUE 

R3  =  FILON(F, ABSTSO, 0. ,1. ,M, 999, 2) 

C 

1  =  1 

DO  55  U  =  0. ,  1. ,  INCR 
KOLU  =  K0L*U 
F(I)  =  MOODIMAG(KOLU) 

1  =  1  +  1 

55  CONTINUE 

R4  =  FILON(F, ABSTSO, 0. ,1. ,M, 999, 2) 

R  :  calculated  value  of  the  integral 


C 

C 

100 


R  =  (R1  -  SGN*R4)  +  J*(R2  +  SGN*R3) 

RR  =  R1  -  SGN*R4 
RI  =  R2  +  SGN*R3 
MAGN2_R  =  RR*RR  +  RI*RI 

ECHO  =  sigina/(2*L) 

ECHO  =  MAGN2_R  *  KOL/8. 

ECHOMO  =  sigma/wavelength 
ECHOMO  =  ECHO*KOL/PI 

WRITE (7,*)  QI,  ECHOMO,  ECHO 
CONTINUE 
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STOP 

END 


REAL  FUNCTION  FILON  (F,  T,  A,  B,  M,  MAX,  KEY) 

C  ********************************************* 

C  - 

C  THIS  PROGRAM  COMPUTES  [F(x)  *COS  (Tx)  ]dx  IF  ''KEY"  =  1, 
C  AND  tF(x)*sin(Tx) ]dx  WHEN  "KEY"  =  2,  USING  FILON'S 
C  METHOD.  FUNCTION  "F(x)"  IS  ASSUMED  TO  BE  GIVEN  AS  A 

C  TABLE  OF  FUNCTIONAL  VALUES  IN  THE  ARRAY  AT  "M" 

C  EQUIDISTANT  POINTS  FROM  "A"  TO  "B",  "M"  ODD. 

C  "MAX"  IS  THE  DIMENSION  OF  THE  ARRAY  DECLARED  IN  THE 
C  CALLING  PROGRAM  (FROM  [REF.  9]). 

C 

C 

INTEGER  KEY,  M,  N,  I,  MAX 

REAL  F(MAX),  AL,  BE,  GA,  SUM,  TH,  H,  S,  C,  SI,  Cl,  S2 ,  C2 
REAL  FI,  F2,  SU,  SUl,  Al,  A,  B,  T,  TH2,  TH4 
N  =  M  -  1 

H  =  (B  -  A)  /  FLOAT  (N) 

TH  =  T  *  H 
S  =  SIN  (TH) 

C  =  COS  (TH) 

TH2  *  TH  *  TH 

IF  (ABS  (TH)  .LE.  0.1)  THEN 
TH4  =  TH2  *  TH2 

AL  =  2.  /  45.  *  TH  *  TH2  *  (1.-  TH2  /  7.  +  TH4  /  105.) 

BE  =  2.  /  3.  +  2.  /  15.  *  TH2  -  4.  /  105.  *  TH4 

GA=4.  /3.  -2.  /15.  *  TH2  +  1.  /  210.  *  TH4 

ELSE 

AL=1.  /TH+S*C/  TH2  -  2.  *  S  *  S  /  (TH2  *  TH) 
BE*(2.  +2.  *C*C  -4.  *S*C/TH)/  TH2 
GA  =  4.  *  (S  /  TH  -  C)  /  TH2 


END 

IF 

SUM 

=  0.0 

FI  = 

F  (1) 

F2  = 

F  (M) 

SI  = 

SIN 

(A 

* 

T) 

S2  = 

SIN 

(B 

* 

T) 

Cl  = 

COS 

(A 

* 

T) 

C2  = 

COS 

(B 

* 

T) 

Al  = 

A  + 

H 

GO  TO  (1, 

2) 

f 

KEY 

SU  = 

F2  * 

S2 

-  FI  * 

SI 

SUl 

=  -0. 

5  * 

(F2  * 

C2  - 

FI  * 

Cl) 

DO  3 

I  = 

2, 

N 

,  2 

SUM 

=  SUM 

+ 

F 

(I)  * 

COS 

(Al  * 

T) 

Al  = 

Al  + 

H 

SUl 

=  SUl 

+ 

F 

(I+l) 

*  COS  (Al 

*  T) 
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3  A1  =  A1  +  H 
GO  TO  4 

2  IF  (T  .EQ.  0.)  THEN 
FILON  =  0. 

RETURN 
END  IF 

SU  =  FI  *  Cl  -  F2  *  C2 

SUl  =  -0.5  *  (F2  *  S2  -  FI  *  SI) 

DO  5  I  *  2,  N,  2 

SUM  *  SUM  +  F  (I)  *  SIN  (A1  *  T) 

A1  =  A1  +  H 

SUl  =  SUl  +  F  (I+l)  *  SIN  (A1  *  T) 

5  A1  =  A1  +  H 

4  FILON  *  H  *  (AL  *  SU  +  BE  *  SUl  +  GA  *  SUM) 

RETURN 

END 

C 

C  ************************** 

REAL  FUNCTION  MOEVREAL(RO) 

C  ************************** 

c 

REAL  K 

COMPLEX  J,  BJC3,  DJE,  FE,  EV,  FACTOR 
COMMON  K,  H,  W,  THETAO,  THETA 
C 

STHO  =  SIN(THETAO) 

CTHO  =  COS(THETAO) 

STH  »  SIN (THETA) 

CTH  =  COS (THETA) 

WROK  =  W*RO/K 
SW  *  SIN (WROK) 

CW  =  COS  (WROK) 

J  =  CMPLX(0. ,1.) 

A  =  -H*W*SW 
B  =  -H*W*W/K*CW 
BA==  ABS(B) 

C  =  SQRT(1.+  A*A) 

D  =  (STHO)* (STHO)  -  (A*CTH0) * (A*CTH0) 

E  =  B*CTH0/(C*C) 

C3  =  C*C*C 
BJC3  =  BA  +  J*C3 
DJE  =  D  +  J*E 
C 

TSE  =  K*H*CW*(CTH  -  CTHO) 

FACTOR  =  2.*C*CEXP(J*TSE) 

C 

FE  =  BA/(2.*C3)*{1-BA/(4.*BJC3))  +  C/ (2 . *BJC3) *DJE 
&  +J*(1+CTH0/C) 

EV  =  FACTOR* FE 
FEVREAL  *  REAL(EV) 

RETURN 

END 
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Q  *********  iHt*  It*  it  *  it  it  it  *  it  ****  it 

REAL  FUNCTION  MOEVIMAG(RO) 

************************** 

REAL  K 

COMPLEX  J,  BJC3,  DJE,  FE,  EV,  FACTOR 
COMMON  K,  H,  W,  THETAO,  THETA 

STHO  =  SIN (THETAO) 

CTHO  =  COS (THETAO) 

STH  =  SIN (THETA) 

CTH  =  COS (THETA) 

WROK  =  W*RO/K 
SW  »  SIN (WROK) 

CW  =  COS (WROK) 

J  =  CMPLX(0.,1.) 

A  *  -H*W*SW 
B  =  -H*W*W/K*CW 
BA=  ABS(B) 

C  =  SQRT(1.+  A*A) 

D  *  (STHO)* (STHO)  -  (A*CTHO) * (A*CTHO) 

E  =  B*CTHO/(C*C) 

C3  =  C*C*C 
BJC3  =  BA  +  J*C3 
DJE  *  D  +  J*E 

TSE  »  K*H*CW*(CTH  -  CTHO) 

FACTOR  »  2.*C*CEXP(J*TSE) 

FE  =  BA/(2.*C3)*(1-BA/(4,*BJC3))  +  C/ (2 . *BJC3) *DJE 
+J*(1+CTH0/C) 

EV  =  FACTOR* FE 
FEVIMAG  *  AIMAG(EV) 

RETURN 
END 

************************** 

REAL  FUNCTION  MOODREAL(RO) 

************************** 


REAL  K 

COMPLEX  J,  BJC3,  DJE,  FO,  OD,  FACTOR 
COMMON  K,  H,  W,  THETAO,  THETA 

STHO  =  SIN (THETAO) 

CTHO  =  COS (THETAO) 

STH  =  SIN (THETA) 

CTH  =  COS  (THETA) 

WROK  =  W*RO/K 
SW  =  SIN (WROK) 
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CW  =  cos (WROK) 

J  =  CMPLX(0.,1.) 

A  =  -H*W*SW 
B  »  -H*W*W/K*CW 
BA==  ABS(B) 

C  =  SQRT(1.+  A*A) 

D  =  2.*A*STH0*CTH0 
E  *  A*B*STHO/(C*C) 

C3  =  C*C*C 
BJC3  =  BA  +  J*C3 
DJE  =  D  +  J*E 

TSE  =  K*H*CW*(CTH  -  CTHO) 

FACTOR  =  2.*C*CEXP(J*TSE) 

FO  *  C/(2.*BJC3)*DJE  +  J*A*STHO/C 
OD  =  FACTOR* FO 
FODREAL  *  REAL(OD) 

RETURN 

END 

****A***********4tA*****4t4t* 

REAL  FUNCTION  MOODIMAG(RO) 

************************** 

REAL  K 

COMPLEX  J,  BJC3,  DJE,  FO,  OD,  FACTOR 
COMMON  K,  H,  W,  THETAO,  THETA 

STHO  «=  SIN  (THETAO) 

CTHO  =  COS (THETAO) 

STH  =  SIN (THETA) 

CTH  =  COS (THETA) 

WROK  =  W*RO/K 
SW  =  SIN (WROK) 

CW  =  COS (WROK) 

J  =  CMPLX(0.  ,1. ) 

A  =  -H*W*SW 
B  =  -H*W*W/K*CW 
BA=  ABS(B) 

C  =  SQRT(1.+  A*A) 

D  =  2 . *A*STHO*CTHO 
E  *  A*B*CTHO/ (C*C) 

C3  *  C*C*C 
BJC3  =  BA  +  J*C3 
DJE  =  D  +  J*E 

TSE  *  K*H*CW*(CTH  -  CTHO) 

FACTOR  =  2.*C*CEXP(J*TSE) 

FO  *  C/(2.*BJC3)*DJE  +  J*A*STHO/C 
OD  =  FACTOR*FO 
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FODIMAG  =  AIMAG(OD) 

RETURN 

END 


H 


PROGRAM  SIMUSBI 
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THIS  PROGRAM  CALCULATES  THE  BISTATIC  ECHO  WIDTH  OF  A 
SINUSOIDAL  SURFACE  z(x)  =  H*cos ( [2*pi/d] *x) 


INPUT  DATA: 

FREQ  =  FREQUENCY  IN  GHz 

INCID  =  ELECTRIC  FIELD  ANGLE  OF  INCIDENSE  (IN  DEGREES) 
WITH  RESPECT  TO  THE  VERTICAL  (+Z)  AXIS 

H,  d  =  PARAMETERS  OF  THE  SINUSOIDAL  SURFACE 

KOL  =  LIMIT  OF  INTEGRATION  (KOL  =  K*L) ,  WHERE  : 

"K"  IS  THE  PROPAGATION  CONSTANT  AND 

"L”  IS  THE  ON  THE  X-AXIS  SURFACE  PROJECTION 

LENGTH 

STEP  =  OBSERVATION  ANGLE  INCREMENT 

M  =  NUMBER  OF  EQUIDISTANT  POINTS  BETWEEN  THE  TWO 
NORMALIZED  LIMITS  (0  AND  1)  OF  INTEGRATION 
(M  ODD,  M  <  999) . 


OUTPUT  DATA: 


ECHOBI  =  NORMALIZED  (WITH  RESPECT  TO  THE  WAVELENGTH) 
BISTATIC  ECHO  WIDTH  AS  A  FUNCTION  OF  THE 
OBSERVATION  ANGLE  THETA. 

ECHO  =  NORMALIZED  (WITH  RESPECT  TO  THE  SURFACE  LENGTH 
2*L)  BISTATIC  ECHO  WIDTH  AS  A  FUNCTION  OF  THE 
OBSERVATION  ANGLE  THETA. 


REAL  F(999) 

REAL  K, KOL, K0LU,INCR, INCID, MAGN2_R 
COMPLEX  J,  R 

EXTERNAL  FILON,  FREAL,  FIMAG 
COMMON  K,  H,  W,  THETAO,  THETAR 

OPEN  (UNIT=8,  FILE='SINUSBI.DAT' ,  FORM= ' FORMATTED ' , 
STATUS= ’ UNKNOWN • ) 

READ(5,*)  FREQ,  INCID,  H,  d,  KOL,  STEP,  M 

PI  =  4.*ATAN(1.) 

W  =  2.*PI/d 
RAD  =  PI/180. 
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THETAO  =  INCID*RAD 

K  =  20.*PI*FREQ/3. 

INCR  =  l./M 
C 

DO  100  QI  =  -90.,  90.,  STEP 
THETAD  =  QI 
THETAR  =  THETAD*RAD 

TSO  =  K0L*( SIN (THETAR)  -  SIN(THETAO)) 

SGN  =  SIGN (1, TSO) 

ABSTSO  =  ABS(TSO) 

C 

1  =  1 

DO  40  U  =  0. ,  1. ,  INCR 
KOLU  =  K0L*U 
F(I)  =  BIEVREAL(KOUJ) 

1  =  1  +  1 
40  CONTINUE 

R1  =  FILON(F, ABSTSO, 0. ,1. ,M, 999,1) 

C 

1  =  1 

DO  45  U  =  0. ,  1. ,  INCR 
KOLU  =  K0L*U 
F(I)  =  BIEVIMAG(KOLU) 

1  =  1  +  1 

45  CONTINUE 

R2  =  FILON(F, ABSTSO, 0. ,1. ,M, 999,1) 

C 

1  =  1 

DO  50  U  =  0. ,  1. ,  INCR 
KOLU  =  K0L*U 
F(I)  =  BIODREAL(KOLU) 

1  =  1  +  1 

50  CONTINUE 

R3  =  FILON(F,ABSTSO,0. ,1. ,M,999,2) 

C 

1  =  1 

DO  55  U  =  0. ,  1. ,  INCR 
KOLU  =  K0L*U 
F(I)  =  BIODIMAG (KOLU) 

1  =  1  +  1 

55  CONTINUE 

R4  =  FILON(F,ABSTSO,0. ,1. ,M,999,2) 

C  R  :  calculated  value  of  the  integral 

R  =  (R1  -  SGN*R4)  +  J*(R2  +  SGN*R3) 

RR  =  R1  -  SGN*R4 

RI  =  R2  +  SGN*R3 

MAGN2_R  =  RR*RR  +  RI*RI 
C  ECHO  =  sigina/(2*L) 

ECHO  =  MAGN2_R  *  KOL/8. 

C  ECHOBI  =  sigma/wavelength 

ECHOBI  =  ECHO+KOL/PI 
C 
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WRITE (8,*)  QI,  ECHOBI,  ECHO 
100  CONTINUE 
STOP 
END 

C  ************************** 

REAL  FUNCTION  BIEVREAL(RO) 

C  ************************** 

C 

REAL  K 

COMPLEX  J,  BJC3,  DJE,  FE,  EV,  FACTOR 
COMMON  K,  H,  W,  THETAO,  THETA 
C 

STHO  =  SIN(THETAO) 

CTHO  *=  COS(THETAO) 

STH  *  SIN (THETA) 

CTH  =  COS  (THETA) 

WROK  *  W*RO/K 
SW  =  SIN (WROK) 

CW  =  COS (WROK) 

J  =  CMPLX(0. ,1. ) 

A  =  -H*W*SW 
B  =  “H*W*W/K*CW 
BA=  ABS(B) 

C  =  SQRT(1.+  A*A) 

D  =  (STHO)* (STHO)  -  (A*CTH0) * (A*CTH0) 

E  =  B*CTH0/(C*C) 

C3  =  C*C*C 
BJC3  *  BA  +  J*C3 
DJE  *  D  +  J*E 
C 

TSE  =  K*H*CW*(CTH  -  CTHO) 

C 

FACTOR  =  2.*C*CEXP(J*TSE) 

C 

FE  =  BA/(2.*C3) *(1-BA/(4.*BJC3) )  +  C/ (2 . *BJC3) *DJE 
&  +J* ( 1+CTHO/C) 

EV  =  FACTOR* FE 
C 

FEVREAL  =  REAL(EV) 

RETURN 

END 

C 

^  ************************** 

REAL  FUNCTION  BIEVIMAG(RO) 

C  ************************** 

c 

REAL  K 

COMPLEX  J,  BJC3,  DJE,  FE,  EV,  FACTOR 
COMMON  K,  H,  W,  THETAO,  THETA 
C 

STHO  *  SIN (THETAO) 
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CTHO  =  COS(THETAO) 

STH  =  SIN (THETA) 

CTH  =  COS (THETA) 

WROK  =  W*RO/K 
SW  =  SIN (WROK) 

CW  =  COS  (WROK) 

J  =  CMPLX(0. ,1.) 

A  =  -H*W*SW 
B  =  -H*W*W/K*CW 
BA=  ABS(B) 

C  =  SQRT(1.+  A*A) 

D  =  (STHO)*(STHO)  -  (A*CTHO) * (A*CTHO) 

E  =  B*CTHO/(C*C) 

C3  =  C*C*C 
BJC3  =  BA  +  J*C3 
DJE  =  D  +  J*E 
C 

TSE  =  K*H*CW*(CTH  -  CTHO) 

FACTOR  =  2.*C*CEXP(J*TSE) 

C 

FE  =  BA/(2.*C3)*(1-BA/(4.*BJC3) )  +  C/ (2 . *BJC3) *DJE 
&  +J*(1+CTH0/C) 

EV  =  FACTOR*FE 
FEVIMAG  =  AIMAG(EV) 

RETURN 

END 

C 

Q  ************A4t******itir**** 

REAL  FUNCTION  BIODREAL(RO) 

C  ************************** 

C 

REAL  K 

COMPLEX  J,  BJC3,  DJE,  FO,  OD,  FACTOR 
COMMON  K,  H,  W,  THETAO,  THETA 
C 

STHO  =  SIN (THETAO) 

CTHO  =  COS (THETAO) 

STH  =  SIN (THETA) 

CTH  =  COS  (THETA) 

WROK  =  W*RO/K 
SW  =  SIN (WROK) 

CW  =  COS (WROK) 

J  =  CMPLX(0.,1.) 

A  =  -H*W*SW 
B  *  -h*W*W/K*CW 
BA=  ABS(B) 

C  =  SQRT(1.+  A*A) 

D  =  2 . *A*STHO*CTHO 
E  =  -A*B*STHO/(C*C) 

C3  =  C*C*C 
BJC3  =  BA  +  J*C3 
DJE  =  D  +  J*E 


79 


o  o 


c 

TSE  =  K*H*CW*(CTH  -  CTHO) 

FACTOR  =  2.*C*CEXP(J*TSE) 

C 

FO  =  C/(2.*BJC3)*DJE  +  J*A*STHO/C 
OD  =  FACTOR* FO 
FODREAL  =  REAL(OD) 

RETURN 

END 

C  ************************** 

REAL  FUNCTION  BIODIMAG (RO) 
************************** 

REAL  K 

COMPLEX  J,  BJC3,  DJE,  FO,  OD,  FACTOR 
COMMON  K,  H,  W,  THETAO,  THETA 
C 

STHO  =  SIN (THETAO) 

CTHO  =  COS (THETAO) 

STH  =  SIN (THETA) 

CTH  =  COS (THETA) 

WROK  =  W*RO/K 
SW  =  SIN (WROK) 

CW  =  COS (WROK) 

J  =  CMPLX(0.  ,1.) 

A  *  -H*W*SW 
B  *  -H*W*W/K*CW 
BA=  ABS(B) 

C  *  SQRT(1.+  A*A) 

D  *  2.*A*STH0*CTH0 
E  =  A*B*STHO/(C*C) 

C3  =  C*C*C 
BJC3  =  BA  +  J*C3 
DJE  =  D  +  J*E 
C 

TSE  =  K*H*CW*(CTH  -  CTHO) 

FACTOR  =  2.*C*CEXP(J*TSE) 

C 

FO  =  C/(2.*BJC3) *DJE  +  J*A*STHO/C 
OD  =  FACTOR* FO 
FODIMAG  =  AIMAG(OD) 

RETURN 

END 
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